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Abstract We introduce a delay differential equation model which de- 
scribes how fish are harvested 
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In our previous studies we investigated the persistence of equation (A) and 
existence of a periodic solution for this equation. 

Here we study the stability (local and global) of the periodic solutions of 
equation (A). 

Keywords-Fishery, Periodic Environment, Delay Differential Equations, 
Global and Local Stability. 



1 Introduction and Preliminaries 

Consider the following differential equation which is widely used in Fisheries 

PE! 

N=[(3(t,N)-M(t,N)}N, (1) 

where N = N(t) is the population biomass, (3(t,N) is the per-capita fecun- 
dity rate, and M(t, N) is the per-capita fishing mortality rate due to natural 
mortality causes and harvesting. 

In equation (JTJ) let P(t,N) be a Hill's type function (UHEl 

M> ir > = Y^Ejr> ( 2 ) 

where a and K are positive constants, and 7 > is a parameter. 

Traditional Population Ecology is based on the concept that carrying 
capacity does not change over time even though it is known 3\ that the 
values of carrying capacity related to the habitat areas might vary, e.g., 
year-to-year changes in weather affect fish population. 

We assume that in © a = a(t) , K = K(t), and M(t, N) = b(t) are 
continuous positive functions. 

Generally, Fishery models |2] recognize that for real organisms it takes 
time to develop from newborns to reproductively active adults. 

Let in equation (J2J) iV = N(8(t)), where 9{t) is the maturation time delay 
< 9(t) < t. If we take into account that delay, then we have the following 
time-lag model based on equation (JTJ) 
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for 7 > 0, with the initial function and the initial value 



N(t) = <p(t), t < 0, N(0) = iV ( 



o 



(4) 



under the following conditions: 

(al) a(t), bit), K(t) are continuous on [0, oo) functions, b{t) > b > 0, K > 
K(t) > k > 0; 

(a2) 9(f) is a continuous function, 9{t) < t, limsup#(t) = oo; 



(a3) ip : (— oo, 0) — > R is a continuous bounded function, (p(t) > 0, N > 0. 

Definition 1.1 A function N : R ^ R with continuous derivative is 
called a (global) solution of problem (J3J), (J3J), if it satisfies equation (jSJ) for all 
t G [0, oo) and equalities (JH) for t < 0. 

If t is the first point, where the solution N(t) of (jHJ), © vanishes, i.e., 
N(t ) = 0, then we consider the only positive solutions of the problem (JSJ), 
dU) on the interval [0,t )- 

Recently [S] we proved the following results: 

Lemma 1.1 Suppose a(t) > b(t), 



Lemma 1.2 Let a(t),b(t), K(t),6(t) be T-periodic functions, a(t) > b(t). 
If at least one of the following conditions hold: 





(bl) 




(b2) 




then equation has at least one periodic positive solution No(t). 
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In what follows, we use a classical result from the theory of differential 
equations with delay |U Ej. 

Lemma 1.3 Suppose that for linear delay differential equation 

x(t) + r{t)x(h{t)) = (5) 

where < t — h(t) < a, the following conditions hold: 

r(t) > r > 0, (6) 

ft 3 
limsup / r(s)ds < - (7) 

t^oo Jh(t) 2 

Then for every solution x of equation^) we have lim^oo x{t) = 0. 

2 Main Results 

Let us study global stability of the periodic solutions of equation 0. 

Theorem 2.1 Let a(t),b(t), K(t),9(t) be T-periodic functions, satisfying 
conditions of Lemma li.il and one of conditions bl) or b2) of Lemma li.M 
Suppose also that 

t 

a{t) > a > 0, 7 J a(s)ds < 6. (8) 
e(t) 

Then there exists the unique positive periodic solution N (t) of (3) and for 
every positive solution N(t) of (3) we have 

limt^iNit) - N {t)) = 0, 

i.e., the positive periodic solution N (t) is a global attractor for all positive 
solutions of (0). 

Proof. Lemma fOl implies that there exists a positive periodic solution N (t). 
If that solution is an attractor for all positive solutions then it is the unique 
positive periodic solution. 

We set N(t) = exp(x(t)) and rewrite equation (3) in the form 

*(*) = - ^(L — -£>(*)• ( 9 ) 

1 + V K(t) ) 
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Suppose u(t) and v(t) are two different solutions of ©• Denote w(t) — u(t) — 
v(t). To prove the Theorem 12. II it is sufficient to show that lim t _ )00 w(t) = 0. 
It follows 



w(t) = a(t) 

Let 



1 1 



,u(0(t))^7 , f e v{0(t))\l 



K(t) J 1 V «■(*) 



and 



7«^) (g|) 7 

c(t) \T 19 5 



M(t) 



(10) 



Using the mean value theorem, we have for every t 

f(y,t)-f(z,t) = f(c)(y-z), (12) 

where min{?/, z} < c(t) < max{y, z}. 
Clearly, 



Equalities (jTTjl - (f*H2|) imply that equation (JTUj) takes the form 

ii>(t) = -M(t)w(6(t)) } (14) 

where 



and 

min{u(0(t)),u(0(*))} < c(t) < m&x{u(9(t)),v(9(t))}. 
Now we want to check that for equation ()14|) all conditions of Lemma fl. HI 
hold. 

From (JT3j) we have M(t) < \^a{t). Therefore inequality (7) holds. Let us 
check inequality (6). Set N±(t) = e u{t \ N 2 (t) = e v ®, where Ni(t),N 2 (t) are 
two solutions of equation (jSJ), corresponding to the solutions u(t) and v(t) of 
equation (JHJ). Lemma fTTTl implies that 

m) > - — v , K _ / 7X2 > o, 



_j_ ^ max{/3 J y 1 ,/3jy 2 } ~j ' 



where and /3n are defined by Lemma fTTTl Hence inequality (6) holds and 
therefore Theorem 12 .11 is proven. 

Consider now equation (jSJ) with proportional coefficients: 
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-br{t) N(t), 



(15) 



where r(t) > ro > 0. Clearly, if a > b then equation (fTKJ) has the unique 
positive equilibrium 

'a 
.b 



N* 



if if. 



(16) 



Corollary 1. If a > b, ( § - 1 ) IT ^ 1, r(t) > r > 0, and 



7alimsup / r(s)ds < 6, 



(17) 



then the equilibrium N* is a global attractor for all positive solutions of 
equation (fTo]) . 

Let us now compare the global attractivity condition (|17j) with the local 
stability conditions. 



Theorem 2.2 Suppose a > b, r(t) > r > and 



j (a — b)b 



~ b b v f MJ 3 
hmsup / r{s)ds < -. 

a t->oo Jo(t) 2 
T/ien i/ie equilibrium N* of equation A 15)) is locally asymptotically stable. 
Proof. Set x = N — N* and from equation (|15|) we have 

ar(t) 



(18) 
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Denote 



Clearly 



F(u,v) 



ar(t) 



^(0,0) = 



6r(t) 



(x(t) + AT*) 



(19) 



j (a — 6)6 



r(t) 



and Fy(0, 0) = 0. Hence for equation (|15[) the linearized equation has a form 

x(t) = - ^ a ~ h ' h r{t)x{6{t)). (20) 

Lemma fOl and condition (fTKJ) imply that equation (j2SJ) is asymptotically sta- 
ble, therefore the positive equilibrium iV* of equation ()15|) is locally asymp- 
totically stable. 

Compare now Theorem 12. II and Theorem 12.21 We have max{6(a — b)} = 
a j '4. Therefore, if 

cry lim sup / r(s)ds < 6, 

t-^oo J0(t) 

then equation ()15|) has locally asymptotically stable equilibrium iV*. 

The last condition does not depend on b, and is identical to condition (JT7|) 
that guarantees the existence of a global attractor. Therefore in Theorem 12. 21 
we obtained the best possible conditions for global attractivity for equation 
Q. 
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